Integrability of the diffusion pole in the diagrammatic description of noninteracting 

electrons in a random potential 
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We discuss restrictions on the existence of the diffusion pole in the translationally invariant di- 
agrammatic treatment of disordered electron systems. We use the Bethe-Salpeter equations for 
the two-particle vertex in the electron-hole and the electron-electron scattering channels and derive 
for systems with time reversal symmetry a nonlinear integral equation the two-particle irreducible 
vertices from both channels must obey. We use this equation to test the existence of the diffusion 
pole in the two-particle vertex. We find that a singularity of the diffusion pole can exist only if it is 
integrable, that is only in the metallic phase in dimensions d > 2. 
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Introduction. Scattering of free charge carriers on 
impurities and lattice imperfections can lead at low- 
temperatures to a metal-semiconductor transition. There 
are two qualitatively different scenarios how a metal can 
turn insulating due to excessive scatterings on impurities. 
In the first case the metal-insulator transition material- 
izes in substitutional alloys when charge carriers are ex- 
pelled from the Fermi surface and an energy gap develops. 
This transition, called split band, is qualitatively well un- 
derstood and quantitatively well modeled by a mean-field 
solution^ The second type of a metal-insulator transi- 
tion is much more complicated and up to now not com- 
pletely understood. Electrons in a metal with random 
impurities can lose their ability to diffuse on macroscopic 
scales. Such scenario was first suggested by Anderson? 
and is now called Anderson localization transition. 

One of principle obstacles of full understanding of the 
Anderson localization transition is inability to describe 
vanishing of diffusion of electrons analytically even in its 
simplest model version and reconcile results from ana- 
lytic and numerical approaches. Analytic, mostly dia- 
grammatic and field-theoretic approaches in the thermo- 
dynamic limit indicate that the critical behavior at the 
Anderson localization transition fits the one-parameter 
scaling scheme with a single correlation length control- 
ling the long-range fluctuations^ On the other hand, 
an increasing number of numerical studies of the Ander- 
son localization transition in finite volumes suggest that 
instead of homogeneous, translationally invariant param- 
eters one has to take into consideration distributions of 
conductances or local particle densities.— The two dif- 
ferent methodological approaches, analytic and numeri- 
cal, disagree not only on the number of relevant control- 
ling parameters needed to understand Anderson local- 
ization but also on the critical behavior and the values 
of the critical exponents. 8 ' 9 Neither of these approaches 
is, however, absolutely conclusive in delivering ultimate 
answers. 

In case of disagreement of results from two rather well 
established and otherwise reliable methods one has to re- 
visit the assumptions under which either results were de- 



rived and to which restrictions they are subject. One of 
the most important features used in the description of the 
critical behavior of the Anderson localization transition 
is a singular low-energy behavior of the density-density 
correlation function of disordered systems. This singu- 
larity has form of a resolvent of a diffusion equation and 
is called the diffusion pole. The existence of the diffusion 
pole and a connection of the diffusion constant with con- 
ductivity are consequences of conservation laws in ran- 
dom systems^ Conservation laws should be a firm part 
of any reliable theory. We, however, showed recently that 
an asymptotic solution of the Anderson model of nonin- 
teracting electrons in high spatial dimensions does not 
fully obey conservation of probability^ We suggested a 
qualitative explanation for such an unexpected behavior 
but more importantly, we amassed arguments that unre- 
stricted compliance with the conservation law is in ran- 
dom systems in conflict with analyticity of the spectral 
function i 12 i 13 Since discussion about the form of the diffu- 
sion pole is still ongoing ; 14 ' 15 we trace down in this paper 
the origin and set exact restrictions on the form of the 
diffusion pole derived within the translationally invariant 
description of disordered systems in the thermodynamic 
limit. We first thoroughly analyze the assumptions used 
to derive the diffusion pole and then prove an assertion 
about the acceptable form of this singularity without re- 
ferring or resorting to any specific approximation. We 
find that in systems invariant with respect to time inver- 
sion the diffusion pole must be integrable in momentum 
space. 

Definitions and assumptions . We model the sys- 
tem of non-interacting electrons by a lattice gas described 
by an Anderson Hamiltonian^ 
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used to capture the impact of randomness on the elec- 
tronic structure of metallic alloys as well as to under- 
stand vanishing of diffusion in the limit of strong ran- 
domness. The first, homogeneous, part of this Hamilto- 
nian is kinetic energy and is diagonalized in momentum 
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space (Bloch waves). The second sum runs over lattice 
sites and describes a site-diagonal random potential. Val- 
ues Vi at different positions are uncorrelated and follow a 
probability distribution P(Vi). This term is diagonalized 
in the direct space by local Wannier states. The two op- 
erators do not commute, quantum fluctuations become 
important and the full Anderson Hamiltonian cannot be 
easily diagonalized. The only way to keep analytic con- 
trol of the behavior of equilibrium states of the Ander- 
son model is to go directly to the thermodynamic limit. 
Standardly it is approached by applying the ergodic the- 
orem, that is, summation over lattice sites equals the 
configurational averaging. This means that we assume 
self-averaging property for all quantities of interest. This 
need not be, however, always fulfilled as we know from 
studies of Anderson localization. Presently we disregard 
this option from consideration as well as the problem of 
the existence of the thermodynamic limit. 

Ergodicity itself, however, does not simplify the pro- 
cess of averaging over randomness. Another assumption 
must be adopted to master this problem. We assume 
that the thermodynamic limit can be performed inde- 
pendently term by term in the expansion in powers of 
the random potential. It means that we expect that the 
configurationally averaged perturbation expansion in the 
random potential converges for all quantities of interest. 

Thermodynamic limit has an important simplifying 
consequence for macroscopic (averaged) quantities. The 
spectrum of a random Hamiltonian in the thermody- 
namic limit is invariant with respect to lattice transla- 
tions. It means that operators H and T R HT^, where 
Tr, is the operator of translation with a lattice vector R, 
have identical spectrum of eigenvalues with translation- 
ally shifted eigenvectors. A lattice translation by a vector 
R„ applied to the Anderson Hamiltonian from Eq. |T]) 
generates a new one, J2k |k) e (k)(k + Y2i N + n Wi{i + n \ 
having the same distribution of random energies. Un- 
less we break translational symmetry in thermodynamic 
states, we are unable to distinguish translationally shifted 
Hamiltonians. We cannot, however, break translational 
invariance of the thermodynamic states arbitrarily, since 
their symmetry should be in concord with the spatial 
distribution of the eigenstates of the Hamiltonian for the 
given configuration of the random potential. Since we do 
not know this spectrum, we must treat all lattice transla- 
tions of the Hamiltonian as equivalent and instead of one 
Hamiltonian we are able to describe only the whole class 
of equivalent Hamiltonians T R HT R . In this way we can- 
not distinguish directly between extended and localized 
eigenstates of the random potential, since the localized 
states are represented by a class of vectors differing by 
lattice translations. 

The natural basis for translationally invariant quan- 
tities is formed by Bloch waves labeled by quasimo- 
menta. We generically denote k, q fermionic and bosonic 
(transferred) momenta respectively. The fundamental 
building blocks of the translationally invariant descrip- 
tion of disordered electrons are averaged one- and two- 
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particle resolvents G(k, z) and G\^, (z + , z_ ; q) , where 
z + = E+uj+ii] and z_ = E—irj are complex energies with 
E standing for the Fermi energy, u> for the bosonic trans- 
fer frequency (energy), and r\ is a (infinitesimally) small 
damping (convergence) factor. We adopt the electron- 
hole representation for the two-particle Green function 
with k and k' for incoming and outgoing electron mo- 
menta. The bosonic momentum q measures the differ- 
ence between the incoming momenta of the electron and 
the hole. Energies of the electron and the hole z+,Z- 
in systems with noninteracting particles are external pa- 
rameters. 

The averaged one-electron resolvent in disordered sys- 
tems can be represented as in many-body theories via an 
irreducible vertex - the self-energy S(k, z). We can write 
a Dyson equation for it 
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The self-energy S(k, z) stands for the impact of the scat- 
terings of the electron on random impurities. Knowledge 
of the self-energy is then sufficient to determine the en- 
ergy spectrum, spectral density and in general all aspects 
of propagation of single particles in disordered media. 

The two-particle resolvent G^ 2 ' can then be repre- 
sented via a two-particle vertex T defined from an equa- 
tion 



G kk'( z +^-;q) 




= G(k, z+)G(q + k, z-) [5{k - k') 
+ r kk ,(^+,^;q)G(k',z + )G(q + k',z_)] (3) 

where <S> denotes the direct product of operators. The 
two-particle vertex introduces a disorder-induced corre- 
lation into the two-particle propagation. Analogously to 
the self-energy it measures the net impact of scatterings 
on impurities on the motion of particles in the presence 
of other particles. 

The two-particle vertex T can further be simplified 
by introducing an irreducible vertex A playing the role 
of a two-particle self-energy. The irreducible and the full 
vertex are connected by a Bethe-Salpeter equation. Un- 
like the one-particle irreducibility, the two-particle irre- 
ducibility is ambiguous. 17 There are two types of two- 
particle irreducibility in systems with elastic scatterings 
only, electron-hole and electron-electron. They are char- 
acterized by different Bethe-Salpeter equations. The 
Bethe-Salpeter equation in the electron-hole scattering 
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channel then reads 



IW(q)=A&(q) 

+ jj E <" (q) G + ( k ")G- (q + k")r 



k"k'(qj 



(4a) 



We suppressed the frequency variables in Eq. (|4a|). since 
they are not dynamical ones. They can be easily deduced 
from the one-electron propagators G±(k) = G R ' A (k) = 
G(k, z±) used there. 

We can introduce another nonequivalent representa- 
tion of the two-particle vertex. If we sum explicitly mul- 
tiple scatterings of two electrons (holes) we can construct 
an alternative Bethe-Salpeter equation 17 

r kk , (q) = A k e k , (q) + iE A kk- (q + k' - k") 



x G+ (k" )G_ (q + k + k' - k")r k » k - (q + k k" 
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We introduced an irreducible vertex in the electron- 
electron scattering channel A ee . Irreducible vertices A eh 
and A ee do not include isolated pair electron-hole and 
electron-electron scatterings, respectively. 

Diffusion pole and electron-hole symmetry. 
Noninteracting particles scattered on impurities are 
marked by a diffusion pole. The low-energy limit 
of a special matrix element of the two-particle resol- 
vent, electron-hole correlation function, has the following 
asymptotics for q — > and u/q — » 
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where np is the density of one-particle states at the Fermi 
levels We used an abbreviation for the energy argu- 
ments Gg£ (E + uj,E;q) = g££ (E + uo + i0+, E - i0+ ; q) . 
The low-energy electron-hole correlation function be- 
comes a propagator of a diffusion equation. 

Such a low-energy behavior is not evident and to prove 
it one has to use Ward identities connecting one- and 
two-particle averaged functions. Ward identities reflect 
conservation laws. In disordered noninteracting systems 
we have probability (mass or charge) conservation. It is 
mathematically equivalent to completeness of the Hilbert 
space of Bloch waves. First Ward identity due to charge 
conservation was derived for disordered systems within 
the mean-field approximation by Velicky 16 and later ex- 
tended beyond this approximation in Ref. IvA It is a 
consequence of an operator identity 
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H z--H 



H 



z--H 



where the multiplication is the standard operator (ma- 
trix) one. This identity holds for any one-particle Hamil- 
tonian. In the thermodynamic limit we must, however, 



average this identity and the averaging procedure need 
not conserve all its aspects when projected onto transla- 
tionally invariant states^ When using the above identity 
in the evaluation of the homogeneous part of the electron- 
hole correlation function, that is q — 0, we obtain 



g,RA 



(0,w) 



-ilj 
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No spatial fluctuations (q ^ 0) of the correlation func- 
tion in the low-frequency limit can be deduced from the 
Velicky- Ward identity. To derive the spatial behavior 
of the diffusion pole in Eq. ([5]) one has to resort to an- 
other relation introduced by Vollhardt and Woffled It 
utilizes the Dyson and Bcthc-Salpctcr equations, Eq. ^ 
and Eq. (|4aj) . and relates the one- and two-particle irre- 
ducible functions S and A eh , respectively. It reads 

£*(q+k,£+c)-I] A (k,£) = i^A^+c^q) 

k' 

x [G R (q + k',E + uj) -G A (k',E)] . (8) 

and was proved diagrammatically (perturbatively) . Us- 
ing the Bethe-Salpeter equation one can show that in the 
homogeneous limit q = this identity reflects the con- 
tinuity equation and hence is equivalent to the Velicky- 
Ward identity. Equation (jHJl together with the Bethe- 
Salpeter equation are then used to show that the long- 
distance fluctuations of the low-energy limit of the cor- 
relation function are controlled by a diffusion constant. 
Introducing a dynamical diffusion constant D(u>) we can 
represent the full leading low-energy asymptotics of the 
electron- hole correlation function as in Eq. ([5]). 1 - Note 
that identity ([7]) holds for both pure and random sys- 
tems, the actual diffusion pole, however, is only the sin- 
gularity from Eq. (J5]) with the momentum dependence of 
the low-energy behavior. To prove such a spatially dif- 
fusive behavior the Bethe-Salpeter equation becomes an 
indispensable tool. 

Another important feature of noninteracting electrons 
on a bipartite lattice without external magnetic field and 
spin-orbit coupling is the time reversal symmetry. Time 
inversion is equivalent to reversing the direction of the 
particle propagation, that is k — > k. The electron 
and the hole interchange their roles. The time-reversal 
invariance for the one-particle propagator then means 
G(k, z) = G(— k, z). Time inversion leads to nontriv- 
ial symmetries when applied onto one of the fermion 
propagators in two-particle functions. The electron- 
hole transformation can be represented either by revers- 
ing the electron line leading to a transformation k — * 
— k',k' — » — k, q — * Q or by reversing the hole propa- 
gator k — > k, k' —¥ k',q — > Q for the electron-hole 
function. Here we denoted Q = q + k + k'. We then 
obtain two symmetry relations for the full two-particle 
vertex 



rkk'(q) = r_ k ,_ k (Q) = rw(-Q). 



(9a) 
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The two-particle irreducible vertices are not invariant 
with respect to time inversion, since the electron-hole 
vertex is transformed onto the electron-electron one and 
vice versa. We then have the following electron-hole sym- 
metry relations 



(9b) 



A kk ,(q)=Af„ h k ,_ k (Q)=A kk ,(-Q). 



This relation says that Bethe-Salpeter equation (|4a|l 
transforms upon time inversion in one particle line onto 
Bethe-Salpeter equation (|4b[) . When the invariance with 
respect to the electron-hole transformation is applied to 
the correlation function we obtain 

(q, u) = ]T G$$ (E + u, E- q - k - k') . (10) 

kk' 

This representation together with the Ward identity, 
Eq. 1(5)), tell us that the same low-energy singularity for 
u, q 2 — » must emerge with the same weight in the aver- 
aged two-particle resolvent G kk (i? + ui, E; q) also in the 
limit w, (k + k' + q) 2 -> 0. 

The uncorrelated propagation of electrons in a ran- 
dom potential does not contain the diffusion pole, and 
hence it must emerge in the vertex function T. Taking 
into account the time-reversal invariance we can single 
out the singular parts of the electron-hole symmetric two- 
particle vertex and obtain 



r^(q,c)= 7k -(q,a; 
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D(uj)q 2 



-iu; + D(uj)(q + k + k') 2 
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The reduced vertex j has a marginal and thermody- 
namically irrelevant singularity for u> — > at k = k' = 
q = 0. It can, nevertheless, display another singular 
behavior in fermionic variables k, k' that is not deriv- 
able from the diffusion pole. Such a singularity must 
not, however, affect the form of the diffusion pole in the 
electron-hole correlation function for q — * 0. The second 
term on the right-hand side of Eq. (fTTj) dominates in the 
leading order of the limit q — > 0, u — > while the third 
one in the limit q + k + k' — > 0, u> — ► 0. We used the dy- 
namical form of the diffusion constant D(oj) so that the 
localization phase would fit. Equation (jlip is the most 
general form of the two-particle vertex reproducing the 
diffusion pole in the correlation function $. The singu- 
larity for q — » is the diffusion pole while the other for 
q + k + k' — » is the Cooper pole caused by multiple 
electron-electron scatterings. To conform this represen- 
tation with Eq. (J7]) we have to satisfy a normalization 
condition that in the metallic phase (-D(O) > 0) reads 



(k)| 



Vkk' 



G+(k')| 2 = 27rn F . (12) 



kk' 



Parquet equations with time-reversal symme- 
try. The full two-particle vertex symmetric with respect 



to the electron-hole transformation can alternatively be 
decomposed by means of the so-called parquet equation 
that can be represented in various equivalent waysi£ 



r kk ,(q) = A kk ,(q) + /C kk ,(q) = A&(q) + /C&(q) 
^,(q)+/C kk ,(q)+J kk (q) 



^kk' 



A k e k '(q)+A&(q)-2kk<(q) 



(13) 



where £ kk /(q) and /C kk /(q) are two-particle reducible 
vertices in the electron-hole and electron-electron chan- 
nels, respectively. We denoted I — A eh n A ee a two- 
particle fully irreducible vertex, that is, a vertex irre- 
ducible simultaneously for both the electron-hole and the 
electron-electron pair propagation (multiple scatterings) . 

The parquet equations hold for the systems where the 
electron-hole and the electron-electron multiple scatter- 
ings are nonequivalent, that is, the corresponding two- 
particle irreducibilities are unambiguous and excluding 
definitions of diagrammatic contributions. The con- 
cept of the parquet theory based on nonequivalence of 
two-particle irreducibility can at best be understood 
in terms of sets of diagrams where addition of func- 
tions is represented by union of sets of diagrams the 
functions stand for. Nonequivalence of the electron- 
hole and the electron-electron multiple scatterings means 
JC ee (~1 IC eh = 0. We trivially have in each a-channel 
A Q n JC a = 0. Further on, we have A eh = A eh n T = 
(A eft n A ee ) U (A eh n /C ee ) ClU/C ee . On the other hand, 

K ee = /c ee n r = {K ee n A eh ) u (/c ee n K eh ) = K ee n A eh . 

Hence K. ee C A eh . Combining the above two relations 
we obtain A eh — I U JC ee from which we reach the par- 
quet representations via irreducible or reducible vertices 
in Eq. ([T3|) , A eh U A ee \I = IU K ee U K. eh = T. 

One must be careful when using the parquet decom- 
position for nonintcracting electrons with elastic scat- 
terings only. In this case multiple scatterings on a sin- 
gle site are identical for both channels. Hence, the two 
Bethe-Salpeter equations (|U) are identical, when the one- 
electron propagators are purely local. We then obtain 
A eh = A ee = T. It means that irreducible and reducible 
local diagrams coincide and the concept of two-particle 
irreducibility becomes ambiguous. To amend this prob- 
lem we introduce a stronger full two-particle irreducibil- 
ity including also local scatterings where the electron and 
the hole are indistinguishable. We denote this vertex J . 
The irreducible vertices X, A eh and A ee for noninteracting 
electrons are then transformed in parquet equations IT3)) 
to 

Xkk'(q) = Jkk'(q) + = J0 ^- J° , (i4a) 



A£ k '(q)=A kk ,(q) + 



1 - J°G+G- 

J°G+G- 
1 - J°G+G- 



J° 



(14b) 



where J° = N~ 3 Ekk'q ^kk' (q) and G± = 
A r_1 ^ k G±(k) are the appropriate local (momentum- 
independent) parts. Vertex A kk (q) is irreducible in 
channel a but does not contain multiple scatterings on 
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the same site. It is important that the fully irreducible 
vertex .Tick' (q) contains only cumulant averaged powers 
of the random potential on the same lattice site so that 
double counting is avoided. 

We now use the symmetries from Eq. ([9]) to replace 
the two irreducible vertices by a single function. We de- 
fine 



A kk ,(q) = A k e k ,(q) = Aft,(-q - k - k' 



(15) 



We use this definition in parquet equation p3[) where we 
represent the full vertex by Bethe-Salpeter equation Pa|) . 
We then obtain a fundamental equation for the irre- 
ducible vertex 

A k k' (q) = 2Tkk' (q) 
+ Jf E Akk " (-q - k - k")G+ (k")G_ (q + k" ) 

k" 

x [A k .. k ,(q)+A k .. k ,(-q-k" - k') - T k -. k ,(q)] . (16) 

This is a nonlinear integral equation for vertex A from 
an input I that may have multiple solutions. We choose 
the physical one by matching it to a perturbative solu- 
tion reached by an iterative procedure with an auxiliary 
coupling constant A and a starting condition A^ ) = AX. 
The iteration procedure for a fixed coupling constant A 
is determined by a recursion formula 



N ^ 



<W - A k r k 7 1} (-q - k - k" )G+(k" ) 
xG_(q + k") (A^ k ,( q )-AJ k .. k ,(q) 

= ^E A k " k " 1) (-q-k-k" ) G + (k") 

k" 

xG_(q|k")A[! k , 1) (-q-k"-k') 



(17) 



In this way vertex A = A^ 00 ) is completely determined 
from the input, the fully irreducible vertex XL. A phys- 
ical solution for A = 1 is reached only if the iteration 
procedure converges for < A ^ 1 and the result can 
analytically be continued to A = 1. This construction of 
the physical solution corresponds to the linked-cluster 
expansion from many-particle physics?^ The iteration 
scheme from Eq. (|17l) is the only available way to reach 
a physical solution and hence its convergence and ana- 
lyticity are of principal importance for the diagrammatic 
description of disordered systems. Using Eqs. (fl4|) we can 
rewrite the above equation to another one for the irre- 
ducible vertex A determined from J . The latter vertex is 
the genuine independent input. Notice that in single-site 
theories with local one-electron propagators we obtain a 
solution A = J° to Eq.([T7|). 

Equation (fl6]) (alternatively Eq. (jTTJ) ) is a fundamen- 
tal equation of motion for the two-particle irreducible 
vertex being electron-hole symmetric. The correspond- 
ing full two-particle vertex obeys simultaneously both 



the Bethe-Salpeter equations in the electron-hole and the 
electron-electron channels, Eqs. ([4]) and the two equa- 
tions are not identical, that is A kk , (q) ^ A kk , (q) . Non- 
linearity of the fundamental equation for the irreducible 
vertex poses restrictions on the admissible form of the 
singular behavior in its solutions. Singularities in the full 
vertex T emerge only via singularities in the irreducible 
vertex A. 

Assertion. Two-particle vertex T of noninteracting 
electrons in a random potential can be decomposed into 
irreducible vertices as 

Ikk' (q) = A kk , (q) + A kk , (-q - k - k') - 2" kk , (q) , 

(18) 

if electrons and holes are distinguishable (non-equivalent) 
quasiparticles and the system is invariant with respect to 
time inversion (electron-hole symmetric). We denoted I 
the two-particle fully irreducible vertex. Irreducible ver- 
tex A obeys Eq. (|16[) . The diffusion pole in the full two- 
particle vertex T may materialize only if it appears in 
the irreducible vertex A. Consequently, the diffusion and 
Cooper poles from Eq. (|11[) can exist in T only in the 
metallic phase in spatial dimensions d > 2. 

Proof. Equation (fT5|) is a direct consequence of 
parquet equation (j 1 3|) where the electron-hole symme- 
try, Eq. (O, is used. The parquet equation holds if the 
electron and the hole arc distinguishable quasiparticles 
via their multiple mutual scatterings. That is, electron- 
electron and electron-hole scatterings do not lead to iden- 
tical results. 

We need not find the most general form of low-energy 
(u) — > 0) singularities compliant with Eq. (fl7j|) but rather 
check whether and when singularities from representa- 
tion (fTTj) can emerge in solutions of Eq. (|16p . 

Vertex A kk < (q) contains the diffusion pole of the full 
vertex T, A kk ,(q) the Cooper pole and the fully irre- 
ducible vertex 2ick'(q) is free of these poles. This conclu- 
sion follows from an alternative form of Eq. (IK 



A kk ' (q) = X kk ' (q) 

+ Jj E Fkk " fa) G + ( k " ) G - fa + k " ) Fk " k ' fa) 



^[A kk "(q)-Z kk "(q)] 



N 



xG+(k")GL(q + k")IW(q) (19) 



where we used the fundamental parquet equation (I13p to 
represent the integral kernel A. The electron-hole sym- 
metry leads in the limit q — > and u> — > to an equation 
for the complex conjugate of the full two-particle vertex 



rfi5(q,a;)*=r^ q)k+cl (-q 3 -a;) 



(20) 



that we use to evaluate the convolution of the diffusion 
poles from the full vertex T in the first sum on the right- 
hand side of Eq. (TT5]). We obtain for k = k' in the leading 
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order of q — -> and 



1 £ If k 4 (q, c)G + (k» )G_ (q + k" )r^ k (q, W ) | 



1 E 



<4G+(k")| 



g^o.u^o TV ^ w 2 + D{u>) 2 q 4 

This squared diffusion pole must be compensated by the 
second sum on the right-hand side of Eq. (fH))) . It means 
that the diffusion pole must be completely contained in 
function Akk' (q) ~2"kk' (q) = ^-kk' (q) ■ From the electron- 
hole symmetry wc then obtain that the Cooper pole must 
completely be contained in function ICfjL, (q) and conse- 
quently the sum of the diffusion and the Cooper poles 
from the full vertex I\k' (q) in Eq- (fTT|) is already part of 
function Tkk' (q) — 2kk' (q) • The fully irreducible vertex 
2kk'(q) is hence free of the diffusion and Cooper poles. 

We discuss first the behavior of the diffusion pole in 
the metallic phase with D(0) = D > 0. When inserting 
the singular part of the two-particle vertex due to the 
diffusion pole we obtain the leading singularity on the 
left-hand side of Eq. (fT7)) 



<Skk'(q,w) 



-y 



N 



—iuj + Dq 2 

^■^G + (k")G-(g + k») 
-?:w + £>(q + k + k") 2 



(21a) 



and on its right-hand side 



SkUq,w) 



-y- 



N jr? -iw + L>(q + k + k") 2 
G+(k")G_(q + k") 



-iuj + D(q + k' + k") 2 



(21b) 



Since the singular term from Eq. (|21aj) contains the com- 
plete form of the diffusion pole, the sum over momenta 
must not bring any new singular contribution in small 
frequencies and is of order O(oj°). To assess the low- 
frequency behavior (w — > 0) of the sum over momenta 
we equal external fermionic momenta k' = k and use an 
asymptotic representation for the contribution from the 
singular part of the integrands 



Skk(q,^) 



^ A q -k^-k, k G J? (q + k)G^(k) 
-iw + Dq 2 
1 A 1 



N ^-f —iw + D(q + k + k" ) 2 ' 
S& (q, w) =p£i q _ k ^_ k , k G*(q + k)G^(k) 



-q- 
1 



Nj-f [-iw + D(q + k + k") 2 Y 



(22a) 



(22b) 



where n is an appropriate momentum cut-off. The 
two expressions cannot be more divergent in the low- 



frequency limit as (— iw)^ 1 for any value of the exter- 
nal momenta q and k. Due to the normalization condi- 
tion, Eq. (TT2"|) , we find to each vector k a set (of mea- 
sure one) of momenta q so that -q-k ^ 0- ^ the 
homogeneous case, q — 0, falls into this set then from 
Eq. (|22ap we obtain integrability of the diffusion pole. 
If not, then for </?k,- q -k 7^ we obtain S kk (q,uj) oc 
{-lujy/^/i-iuo + Dq 2 ) and Si\{q,w) oc {-iw) d l 2 ~ 2 . 
For low dimensions d < 2, both functions S kk (q,u>) and 
S kk (q,uj) have a stronger divergence than (— iuj)^ 1 (for 
q = 0) and Eq. (fTT|) cannot be satisfied by any func- 
tion Akk(q, w). The diffusion pole can hence exist in the 
metallic phase only in dimensions d > 2. 

In the localized phase we expect the following low- 
energy asymptotics (q — > 0, u> — > 0) of the dynamical 
diffusion constant 



= — — > 



(23) 



where £ is a localization length. Using this asymptotics 
we can represent the singular part of the irreducible ver- 
tex A as follows 



A sinn i \ 
kk' (q> w ) 



<^kk' 



i 



-iw I + £ 2 q 2 



(24) 



We utilize the electron-hole symmetry to evaluate the 
complex conjugate of the irreducible vertex A in the low- 
frequency oj — > and momentum q — > limit 



\ RA 



k'-ql 



-w) 



(25) 



and use it to derive a condition for vanishing of quadratic 
singularity of order ui~ 2 on the right-hand side of 
Eq. (fTB|) . After substituting the representation of ver- 
tex A from Eq. (|24]) and setting k' = k and q = in 
Eq. (|16[) we obtain 



N 4^ 



<k-G + (k") 



1 +£ 2 (k + k") 2 



[2 + C 2 (k + k") 2 ] =0. (26) 



This condition can be fulfilled only if the irreducible ver- 
tex Akk" (q, w) is free of the singularity due to the dif- 



fusion pole for q — > and w — > 0, that is, <p kk , 
point-wise. The diffusion pole hence cannot exist in the 
localized phase. 

Discussion and conclusions. The most severe 
consequence of the Assertion is nonexistence of the dif- 
fusion pole in the localized phase in any dimension. It 
means that when approaching the low-energy limit q — > 
and io I q — ► in low dimensions [d < 2) we cannot meet 
the diffusion-pole singularity. The localized phase must 
be reached in a non-critical or a less critical manner than 
that of the diffusion pole. Theories, such as the self- 
consistent theory of Anderson localization of Vollhardt 
and Wolfle, 20 leading to solutions with a nonintegrable 
diffusion pole are in conflict with the Bethe-Salpeter 
equation either in the electron-hole or in the electron- 
electron channel or with the electron-hole symmetry at 
the two-particle level. 
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The Assertion poses no restriction on the expected 
form of the diffusion pole in the metallic phase in di- 
mensions d > 2, since the singularity is integrable. The 
localized phase in d > 2 is, however, different. There 
the widely accepted behavior of the diffusion pole, due 
to vanishing of the diffusion constant (D = 0), becomes 
momentum independent and hence nonintegrable. The 
fundamental equation (|f 6[) for the irreducible vertex A 
cannot lead to a two-particle vertex with such a singular- 
ity. If the diffusion pole in d > 2 survives in the metallic 
phase unchanged till the Anderson metal-insulator tran- 
sition, there must be a jump at the transition point at 
which the diffusion pole abruptly ceases to exist. 

Numerical simulations nevertheless seem to confirm 
the existence of the diffusion pole in the localized phased 
There are two possible conclusions we can draw from 
these incommensurable results. One can speculate that 
some of the assumptions on which the diagrammatic 
translationally invariant description of random systems is 
based do not hold near to the Anderson localization tran- 
sition. Either ergodicity may be broken or one cannot 
average the expansion in the random potential term by 
term, or the concept of distinguishability of the electron- 
electron and electron-hole scatterings is invalid. If this 
was true, then one had either to revisit the derivation of 
the diffusion pole, being presently heavily based on the 
Bethe-Salpeter representation of the two-particle vertex, 
or to question the concept of electrons and holes as dis- 
tinguishable quasiparticles in the localized phase. 

On the other hand, numerical simulations are per- 
formable only on rather small lattices where one can- 
not effectively reach the diffusive regime q — > with 
uj/q — > 0. Ref. HI investigates the opposite limit w — > 
with q/ui — ► 0. As we know^ the two limits do not com- 
mute and the latter has no relevance for the existence of 
the diffusion pole. The numerically observed 1/lu behav- 
ior reflects only the Velicky identity valid for random 
as well as pure systems. A conflict between the form of 
the diffusion pole and the Bethe-Salpeter equation arises 
only in the critical region of the latter. Ward identity |7]) 
cannot be extended to inhomogeneous long-range fluctu- 
ations and the homogeneous low-frequency limit to — > 
with q = may be a singular point having no macro- 
scopic relevance for nonzero spatial fluctuations q > 
in the thermodynamic limit. 

Last but not least, we obtain as a consequence of 
Eq. (jf 6p that the two-particle vertex rkk'(q) in the 



metallic phase of the most interesting spatial dimen- 
sions 2 < d < 4 contains apart from the diffusion and 
the Cooper pole also another low-energy singularity for 
lu -> and |k - k'| -> 0. We found that S^(0,w) = 
(— iLo) d / 2 ~ 2 in d < 4 and hence a new singularity in 
vertex A k ? j^(q, uj) emerges for k — k' — > 0. Due to the 
normalization condition, Eq. (jf 2[) . it must be integrable, 
which is the case for d > 2. This new singularity is com- 
patible with the decomposition from Eq. (jllj) of the two- 
particle vertex T into singularities caused by the diffusion 
pole. The existence of a new singularity makes either the 
weight of the diffusion pole (p^, or 7 k ^'(q, cj) or both 
singular with an integrable singularity. A new singular- 
ity in the two-particle vertex indicates that the averaged 
two-particle functions in spatial dimensions 2 < d < 4 
behave qualitatively differently and have a richer ana- 
lytic structure from those in higher dimensions. How far 
this singularity influences the macroscopic behavior and 
transport properties of disordered systems and in par- 
ticular criticality of the Anderson localization transition 
remains to be investigated. 

To conclude, we proved in this paper that the diffusion 
pole in the two-particle vertex can exist in the models of 
nonintcracting electrons in a random potential with time 
reversal symmetry only in the metallic phase in dimen- 
sions d > 2. An equation of motion for the two-particle 
irreducible vertices prevents the existence of the diffusion 
pole in the localized phase. The existing translationally 
invariant descriptions of electrons in a random potential 
predicting the existence of a pole in the localized phase 
should hence be revisited. In view of our result, it seems 
very difficult, if not impossible, to build up a consistent 
analytic theory of Anderson localization with the diffu- 
sion pole in the localized phase. 
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